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DIGITAL COMPUTER 



LIBRARY KOUTIKE L 6 - 22k 



TITLE 

TYPE 

NUMB3HS OF WORDS 

TEMPORARY STORAGE 



Solution of a Set of n Simultaneous Linear Algebraic 



Equations using Magnetic Drum Storage 

Open 

162 

9 words at thru 8 

n+1 words at r~, r, , .»., r 

n+2 words at w Q , w., . .., 

^i~2i words at k^, k^, 



n* n+1 



• • • 



11* 



on J 



In' 



(DOI or SADOI) 



CAPACITY 
ACCURACY 
PARAMETERS 



k , s . 
n-l,n, n-1 

Where k (i < j ) is the location of a ,, an element of 

the triangularized matrix 

and s. is the location of the sum check of the (i+l^st row 

of the triangularized matrix 

1^3 equations 

Depends on the conditioning of the equations . 

S3* Sk y . . ., SS. During input of the routine, the following 

parameters are required 



3: 


OOF 


00 r F 


k: 


OOF 


00 r F 
n 


5: 


OOF 


00 nF 


6: 


OOF 


OOk^F 


7: 


OOF 


00 iF 


8: 


OOF 


00 v F 


9: 


OOF 


00 v F 
n 


K? 


OOF 


00 eF 


S: 


OOF 


00 dF 


where 


r o' 


r -j t • • 



. r are used for the results x~, x. . 
n 1 



. . ., x (see description); 



n 



k^Q is the first location of the reduced triangular matrix 
in the electrostatic memory^ 
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DESCRIPTIONS 



SCALING THE DATA; 



i is the address of the first order of the input subroutine; 

w 0' w l* * * •* w n ' w n +l are used for "t^Porary storage; 
e is the location to which control is transferred when the 
routine is finished, 
d is the first location of the drum to be used. 

This routine will read a set of equations, reduce them to 
triangular form, and then calculate the solution vector 
x 0* x l> • ' •* x n _i* If tne triangular matrix exceeds 
location 989, then the successive elements are stored on 

the drum beginning at d. * is placed in location r = r +n. 

n n 

x n is the scale factor (a power of ten) such that the roots 

are x Q /x n , X-^x^ . . ., x^-j/x^ Although the scale factor 

will be a power of l/lO one may alter this property of the 

routine. If the word in U7, relative, is made l/r instead 

of l/lO the scale factor will be power of l/r. 

To solve the equations 



n-1 



j=0 



x. + a. 
in 



■0, i = 0, 1, 



o . . 



, n-1 



each of the elements a 4 must be scaled to be less than one- 
half in absolute value. Scaling of the a. . matrix by rows and 
columns should be done so as to make the magnitudes of the 
elements as large as possible consistent with having them 
less than 1/2. 

DISCUSSION OF THE METHOD The method used here is the one used in library 

routine LI. Several slight changes have been made to keep the 
maximum amount of precision. First, when an element of a rev 
exceeds one-half, the elements are all multiplied by 4/10 
rather than by l/lO as in LI. Secondly, before a row becomes 
part of the triangular matrix, the row is scaled up so that 
the largest element in the row will become greater or equal 
to a fourth but less than a half. 

SUM CHECKS 2 If a row has been read from the drum inaccurately during the 
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triangularization, then the program will stop on an FF 
instruction. The right hand address will be (26) /-. If 
the sum check fails during back substitution, the right 
hand address of the FF instruction will be (27). g. The 
back substitution sequence can be begun again by flipping 
the white switch. 



lgr 



DA3 B Oct»by 2k, 19g6RT; 5/20/5 $ 
CODED B Y y^ « AdLL- ~ 
APPROVE) W i^)f t ?#» 



I r 

! 




...JQRDBBL.-.- 
00 K(L6)i 
i+l 3F 




TiQXES. PAGE 1 








hi 6f 




Clear Counters 


l 


in 8F 








U k*F 


From 86 


Set i = 


2 


50 S3 








50 2L 




Input row u, (a^, V ' ' *' a^) into 


3 


26 S7 




(r Q , r 1 ,...,r n-1 ) 




L5 150L 




m = 0, 1, . • » , n-1 


k 


L5 133L 






5 


1+0 2.1L 








1*1 7F 


from 82,84 




6 


L5 135L 




Set initial values 




iA 1*F 




w = w Q 


7 


1*2 73L 
U6 33L 




(for triangularization) 

i 


8 


1*2 22L 
1*2 35L 




1 


9 


1*6 3% 
1*6 i*3L 




i 


10 


1*2 i*7L 
1*2 67L 






11 


1*2 69L 
k6 19L 






12 


1*2 28L 
1*9 2F 






13 


L5 136L 
Lh 1*F 




Set initial values r = r Q 


ll+ 


1*0 58L 

1*2 33L 




(for triangularization) 


15 


k6 1*5L 
k6 37L 






16 


1*2 3^L 

1*2 39L 

. ... 







l6 



LOCATION 



17 
18 

19 
20 
21 
22 

23 
2k 

25 
26 

27 
28 

29 
30 
31 
32 
33 
3k 



ORDER 



k2 Jul 

L5 13^L 
42 kjL 
k6 k2h 
4l (w)F 
L3 8F 
50 1023S4 
36 33L 
JO F 

L5 (a Q0 )F 
kO. F 
40 (w)F 
F5 21L 
40 21L 
LO 137L 
36 26L 
L5 138L 
40 21L 
F5 22L 
40 22L 
LO 152L 
36 21L 
L5 F 
LO (w)F 
kO F 
F5 28L 
kO 28L 
LO 161L 
36 28L 
L3 F 
36 33L 
FF 38F 
L3 (w/F 
L6 (r)F 
36 37L 

L5 (r)F 



NOTES 
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L 6 



by 11 



From 27 
by 5,23,25 

by 8,26 



Set both addresses in interchange 
order to 1 



— If last row, jump to 33L 



Read row of . triangular i zed matrix 



into w. 







w 



1* 



n-1 



from 2k 



from 31 
by 12,30 



by 1, k6 
from 20,32 
by Ik 



by 16 



Test to see if storage will exceed location 989 
in electrostatic memory 



If sum check fails, then stop 



Test |N(r k )| - |N(w k ) 



35 

36 

37 
38 

39 

40 

in 

^3 
44 

45 

46 

47 
48 

49 
50 

51 
52 



QRP5R 
50 6f 
66 (w)f 
47 42L 
22 40L 
50 (r)F 
S3 F 
36 40L 

50 l40L 
75 (w)F 
66 (r)F 
43 43L 

51 F 
40 5F 
L5 (r)F 
40 (l)F 
50 5F 
L5 (w)F 
40 (l)F 
7 J IF 
L4 F 

40 (r)F 
L4 2F 

56 47L 
47 33L 
L5 IF 
40 (w)p 
L3 IF 
L6 7F 
32 50L 
L7 IF 
40 7F 
F5 4lL 
42 4lL 
00 20F 
46 45L 
F5 47L 



NOTES 
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L 6 



by 8 



by 15,925 
from 34 



Divide and arrange 
no interchange 



by 9 
by 16 

from 36 



by 17,51 
from 55 
by 18,36 



by 9,54 
by 18, 40 



by 15,52,93 



from 46' 
by 10,53 



from 49 



Test for zero 

Take care of case when w.=r 

Divide and arrange 



Store ratios = -k in 5 



if new r is < -l/2 or > +l/2, this fact is 
recorded by clearing address of 33 • 



K. -* v ( no interchange ) 

or r. -» w . (interchange) 
1 i x ' 

Determine largest element of row 



Step various addresses 
dependent on i, and loop until 
after i = n has been done. 













r 

LOCATION 


ORDER 


i NOTES PAGE 4 




53 
5* 

55 
56 


42 471 
00 20F 
46 k% 
LO 139L 

32 Ml 

4} 68L 
L5 33L 
LO 48L 








57 


52 61L 
50 l4lL 


from 6l 






58 


7J (r)F 
HO (r)F 


by 14,60 
by 14,60 






59 


L5 58L 
L4 134L 








60 


40 58L 
LO 142L 


• 






61 


3? 57L 
L5 ^2L 


from 57 






6? 


LO 132L 
36 1531 




Determine if interchange has taken place ^ 




63 


50 6F 
F5 68L 


from 66 






61* 


42 68L 
L7 IF 








65 


00 IF 

kO 7F 








66 


32 63L 

41 1S9 








67 


00 17 
L5 (w)F 


by 10, 71 


,73 Scale up elements 




68 


10 2F 
00 ( )F 


by 55,64 


< 




69 


32 69L 




- 




70 


kO (w)F 
Ik 1S9 
kO 139 


by 11,72 


Coapute sun check 



L6 



LOCAJTIOI 


asam 


ICOBS EM»5 


71 


F5 67L 

Ho 67L 


- 


• 


72 


k2 69L 
L0 1^3L 




! 


73 


32 67L 








L5 (v)F 


by 7,78* 


fraa 79 


7* 


/JO F\ 
[kO F/ 


by lt.,75, 77,155*157 Store row if interchange has taken place. 








75 

- 


F5 7^L 
*0 jkL 






76 


LO l^L 
36 78L 






7T 


L5 1^5L 
1*0 7^L 






78 


F5 73L 


froa 76 






k2 73L 






79 


LO 151L 
32 73L 






80 


L5 ^F 
IA 13^L 


fraa 158 




81 


kO kF 
LO 3F 




, 


82 


56 83L 
22 5L 




• 


83 


1*0 8F 
L3 8F 


fraa 82 


• 


8U 


32 5L 

L5 5F 








lA 15^L 
40 3F 




1 


i 86 

i 
i 


LO lk6L 
52 1L 




i 


! 8j 


L5 l^TL 


frcK H5 






kO Sk 


fraa 121 


x^ * 1/10 


88 ' 


19 37F 
1*0 6f 







I LOCATION 



89 

90 

91 

92 

93 

9* 

95 

96 

97 

98 

99 

100 

101 

102 

103 

104 

105 

106 

107 



GRDM 



7k 21L 
40 7F 
L5 20L 
46 127L 
L5 70L 
LO 134L 
LO 134L 
46 37L 
46 126L 
46 4% 
46 124L 
L5 120L 
42 117L 
L5 115L 
46 117L 
L5 7F 
FO 6F 
LO 138L 
32 99L 
F4 l48L 
26 100L 

L4 138L 

40 7F 
4o 10JL 
L5 47L 
LO 6f 
42 110L 
42 104L 
^85 11F\ 
^00 F / 
40 F 
40 (w)F 
F5 105L 
40 103L 
LO 137L 
56 108L 
L5 138L 
40 103L 



HOEEB 



PUB 6 



I* 



from 131 



from 98 
from 99 



from 109 

by 100,105,107 



by 102 



LOCATION 



108 
109 
110 
111 
112 
113 
114 

115 
116 
117 
118 

119 
120 
121 
122 
123 
124 

125 



H- 



F5 1C&L 
42 104L 
LO 152L 
36 103L 
L5 F 
LO (w)F 
40 F 
F5 110L 
42 110L 
LO 161L 
36 110L 
L3 F 
36 ll6L 
FF 39F 
50 S9 
26 87L 
50 F 
S5 F 
50 (w)F 
74 (r)F 
L4 F 
40 F 
L4 2F 
32 121L 
50 147L 
7J S4 
22 87L 
L5 117L 
LO 134L 
40 117L 
LO 37L 
32 116L 
L3 (w i )F 
L6 F 
36 120L 

L5 F 



bokbs 



from 106 



from 113 
by 102,112 



BMS 7 



JLi. 



Test sum check 

If sum check has failed, back substitution 

begin again by flipping white switch. 



can 



from 114,123 



by 96,122 1 Evaluate scaler product 
by 95,122 U v x 
s=i+l 8 S 



from 125 j If we are liable to exceed capacity 



reduce x„ by 10 
n 



from 119 



by 94 



LOCATION ORDER 

126 66 (v t )T *y'93 

81 F 
12J 40 (r t )* *7 90,128 

L5 127L 

128 LO 134L 
46 127L 

129 LO l49L 
56 SK 

150 F5 6f 

40 6F 
131 L5 37L 

26 92L 

132 ho if 

00 F 

133 JO F 
L5 86 

134 00 IF 
00 IF 

135 LTS8 
12 88 

136 V S3 
40 S3 

13T 50 F 

L5 99QF 

138 85 11F 
00 88 

139 MO 189 

■ ; 00 f 

140 7L 4095F 
LL 4095F 

141 00 F 
00 1*000 0000 OOOOJ 

l4fi U 184 

40 184 
143 80 IF 

L5 189 



10TB 



FlfflE 8 



LS 



llak 



"location 



144 

145 
146 
147 
148 
1^9 
150 

151 
152 

153 
154 

155 
156 

157 
158 

« 159 

i 

I 160 

I 161 



II I i n ii piu 'i 



ORDER 



NOTES 
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L6 



50 F 

40 990F 

86 11F 

00 SS 

80 S5 

00 S5 

00 F 

00 1000 0000 OOOOJ 

JO F 

L5 989F 

NO S3 

L5 F 

JO F 

40 S6 

S2 6JL 

L5 2S9 

NO F 

40 2S9 

L5 159L 
LO 4F 
10 20F 
L4 74l 
40 7^L 
LO 144L 
36 I58L 
L4 160L 
L4 145L 
40 74L 
26 80L 
00 F 
00 2S5 
00 2S5 
80 F 
00 F 
75 F 
LO 1S9 



from 62 



f rom 156 



